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MILLIMETER-ACCURACY 
SATELLITE  NAVIGATION 


Definition: 

Determining  instantaneous  position 
coordinates  of  a  point  on  a  moving 
platform  (land  /  sea  /  air  /  space 
vehicle)  with  respect  to  a  reference 
frame  defined  by  points  fixed  on  the 
ground,  with  1-mm  uncertainties. 

“Satellite  navigation”  means  navi¬ 
gation  using  satellites.  (The  moving 
platform  may  also  he  a  satellite.) 
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MILLIMETER-ACCURACY 
SATELLITE  NAVIGATION 

Object  of  this  talk: 

To  convince  you  that  mm-level 
navigation  is  possible. 


Plan: 

1 .  Show  you  how  positions  of  fixed 
points  have  been  determined  with 
millimeter  accuracy. 

2.  Show  you  how  positions  of  moving 
points  have  been  determined  with 
centimeter  accuracy  under  certain 
conditions. 

3.  Discuss  problems  and  solutions. 
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Re:  HOW  POSITIONS  OF  FIXED 
POINTS  HAVE  BEEN  DETERMINED 
WITH  MILLIMETER  ACCURACY 


1.  The  basic  method. 

2.  Some  results. 

3.  What  the  major  problems  have 
been,  and  how  they  have  been 
solved. 
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BASIC  METHOD  (simplified): 


1 .  Observe  LI  -  &  L2-band  radio  signals 
transmitted  by  GPS  satellites. 

2.  Use  observations  from  3  fixed  ref¬ 
erence  points  to  determine  satellite 
orbits  with  respect  to  the  coordinate 
frame  defined  by  these  points. 

3.  Use  observations  at  an  unknown 
position  to  determine  the  unknown 
position  coordinates. 
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TRICK  #1: 


Observe  &  use  the  carrier-wave 
phases  of  the  GPS  signals. 

Easy  to  measure  within  0.1  radian, 
equivalent  to  3  mm  of  path  length. 

Possible  to  measure  within  0.03  rad, 
equivalent  to  1  mm. 
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TRICK  #2: 


Use  only  doubly-differenced 
carrier-phase  observations. 

Difference  between  simultaneous 
observations  of  same  satellite  at  two 
receiving  points  cancels  sensitivity  to 
satellite-transmitter  phase. 

Difference  between  simultaneous  ob¬ 
servations  of  any  two  satellites  at  any 
receiving  point  cancels  sensitivity  to 
receiver-related  phase. 
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Double-Differencing  (continued): 


In  general,  sensitivity  to  any  “common¬ 
mode”  effect,  including  any 

•  receiver-related  (satellite-independent) 
effect,  or  any 

•  satellite-related  (receiver-independent) 
effect, 

cancels  when  observations  are  doubly- 
differenced. 

The  only  surviving  effects  are  those 
which  are  “doubly  different,” 

i.e.,  different  for  observations  of  different 
satellites,  differently  at  different 
receivers. 
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Surviving  Effects,  Bad  &  Good: 


Chief  of  the  baddies  (error  sources)  is 
tropospheric  refractivity.  E.g.  at  one  of 
the  receiving  points  but  not  the  other, 
one  of  the  sate'lites  but  not  the  other 
appears  behind  a  cloud.  A  problem  if 
observing  from  under  the  weather. 

The  desired  sensitivity  of  the  observa¬ 
tions  (to  the  unknown  position)  is  not  at 
all  reduced  by  differencing  between  this 
position  and  another,  and  is  not  much 
reduced  by  differencing  between  satel¬ 
lites  if  these  satellites  are  well  separated 
in  the  sky. 
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FIXED-POSITIONING  DEMONSTRATION 


COMMENTS: 

1 .  All  antennas/receivers*  removed  and 
replaced  each  time. 

2.  All  receivers  were  codeless  (immune 
to  GPS  “A-S”  and  “S-A”). 

3.  Nothing  but  crystal  oscillators. 

4.  No  observations  except  from  the  1 2 
points  shown  (3  ref.  +  9  unk.  pts.)** 

5.  Daytime,  solar  maximum  (high  and 
rapidly  varying  ionospheric  refraction). 

6.  Gulf  coast  weather  (humid,  lots  of  clouds) 


*  except  at  two  of  the  reference  points,  at  Austin  and  Nacogdoches. 
“Three  other  stations  were  used,  on  one  day,  to  determine 
coordinates  of  reference  points. 
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WHY  IT’S  FUNDAMENTALLY  HARDER 
TO  DETERMINE  POSITION  ON  A 
MOVING  PLATFORM  (AS  OPPOSED 
TO  A  FIXED  POINT): 

Errors  can’t  be  time-averaged.  (The 
worst  error  source,  doubly-different 
tropospheric  refraction,  varies  over 
minutes  to  hours.) 


WHY  IT  MAY  (OR  MAY  NOT)  BE  HARDER 
IN  PRACTICE: 

Resolving  the  integer-cycle  ambiguities 
of  the  doubly-differenced  phase  obser¬ 
vations  is  harder  if  you  can’t  wait  for  the 
satellites  to  move  to  different  positions 
in  the  sky. 
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AMBIGUITY  RESOLUTION 


The  interpretation  of  a  (doubly-differenced) 
phase  observation  in  terms  of  position  is 
ambiguous.  The  “likelihood”  function  of 
position,  given  a  set  of  phase-difference 
observations  Act  of  one  or  more  satellites, 
differenced  between  the  unknown  position  r 
and  some  reference  position,  is 


=  y  e2m  (  Atfefc  -  r .  sjkA ) 

over  all  over  all 
observ’n  satellites,  k 
times,  j 


(  Sjk  =  unit-vector  in  direction  of  kth  satellite  at  jth  observation 
time; 

X  =  carrier  wavelength;  sum  may  include  observations  at  multiple 
wavelengths.) 
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AMBIGUITY  EXAMPLE: 


Suppose  that  there  is  just  one  observation 
time  (e.g.  because  we  want  instantaneous 
position  on  a  moving  platform)  and  that  two 
satellites,  in  directions  ^  and  s2 ,  are 
observed  at  this  time.  Then  the  position- 
likelihood  is  maximized  by 

r.(srs2)=  X( Afo  - Afc  +  n  ) 

where  n  is  any  integer.  The  maximum- 
likelihood  position  estimate  is  the  locus  of 
points  on  an  infinite  set  of  parallel  planes, 
perpendicular  to  (s^-s2)  and  separated 
along  this  direction  by  A/|srs2|. 
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ANOTHER  EXAMPLE 


Six  satellites  observed  at  LI  wavelength  only, 
for  15  min,  (each  satellite  moves  ~  1/8th  rad.). 


LIKELIHOOD  FUNCTION  (in  E-W  plane): 


EAST  WEST 

1  ■  8  1  a 


Percent  of  Max. 
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SAME  EXAMPLE,  except  obs’ns  continue  for 
two  hours  (each  satellite  moves  -1  rad.). 


LIKELIHOOD  FUNCTION  (in  E-W  plane) 


Fig.  2.  Like  Fig.  1,  except  the  period  of  observation  is  2  h  instead  of 

15  min. 
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ANOTHER  EXAMPLE  (provided  by  G.  L. 
Mader  of  NGS,  Rockville): 

•  7  satellites 

•  LI  &  L2  wavelengths 

•  1  point  in  time 

•  moving  platform  (P-3  airplane  nr.  Pax  R.) 

LIKELIHOOD  FUNCTION  (3-D,  perspective 
view) : 


P-3  PaxKluer  12/12/89  1  =  In  L1«L2  1  pt  11:34:38  b  =  178 lZw  alrbc 


FOR  INSTANTANEOUS  AMBIGUITY 
RESOLUTION: 

•  observe  both  LI  &  L2 

•  account  for  ionosphere  (unless  very  near 
reference  point) 

•  observe  at  least  7  satellites 

OPTIONS: 

•  resolve  ambiguity  once  and  maintain  track 

•  combine  GLONASS  &  GPS 

•  use  other  satellites  (MARISAT ; 
GEOBEACON) 
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[57]  ABSTRACT 

Techniques  are  disclosed  for  determining  orbital  data  of 
space  borne  vehicles  including  earth  satellites  such  as 
those  of  the  NAVSTAR  Global  Positioning  System 
Each  of  a  set  of  such  satellites  transmits  signals  which 
include  earner  waves  which  may  be  suppressed,  or  only 
implicity  present.  The  signals  are  received  from  the 
observable  satellites  concurrently  by  means  of  an  an¬ 
tenna  at  each  of  at  least  three  ground  stations  forming  a 
network  of  baselines  The  stations  are  arrayed  such  that 
the  ratio  of  the  maximum  lo  (he  minimum  baseline 
length  is  much  greater  than  one.  From  the  signals  re¬ 
ceived  at  a  station  pair  forming  each  baseline  a  time 
senes  of  doubly-difTerenced  phase  measurement  data  is 
formed  which  is  biased  by  an  integer  number  of  cycles 
of  phase.  The  data  senes  for  different  satellite  and  sta¬ 
tion  pairs  are  processed  together  to  determine  the  orbits 
of  the  satellites  and  the  doubly-difTerenced  phase  biases 
Unique  determination  of  the  integer  values  of  at  least 
some  of  the  biases  is  facilitated  by  the  above  noted 
spatial  arrangement  of  the  stations  such  that  the  ratio  of 
the  maximum  to  the  minimum  baseline  length  is  much 
greater  than  one  This  integer  bias  determination  en¬ 
hances  the  accuracy  of  the  related  orbtt  determination 
Unique  determination  of  the  integer  values  of  at  least 
some  of  the  doubly-differenced  earner  phase  biases  may 
also  be  facilitated  by  the  use  of  a  plurality  of  earner 
frequencies  with  the  ratio  of  the  maximum  to  the  mini¬ 
mum  frequency  being  much  greater  than  one. 

18  Claims.  6  Drawing  Sheets 
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